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1. Introduction 

The well known NKG function is a very useful tool to describe 
the lateral extension of the electromagnetic component in EAS ; 
however, in spite of non negligible qualities (simplicity, normali- 
zation by beta-funct ion) , it doesn’t correspond exactly to the 
natural shape of the lateral electron distribution. Several bias may 
occur in size estimation if NKG is used without correction, for 
instance, contradiction between lateral and longitudinal development 
the lateral parameter Sj^ being quite lower than the longitudinal one 
(1)(2). We emphasize here how the longitudinal age parameter s can 
be correlated with the information obtained from the lateral elec- 
tron densities according to the conditions of use of the NKG 
function. 


2. Local age parameter phenomenology 

The theoretical age parameter s^ illustrates the declining stage 
of a shower and is determined as the saddle point in the inverse 
Mellin transformation. We have postulated that the NKG function, 
derived from diffusion equations fails to describe EAS data, mainly 
because a uniqueness of parameter s was assumed, and we have admit- 
ted local agreement with NKG function in small bands of distance by 
introducing the local age parameter 


= 2 ^ X + 2)) if = NKG) 


(x = r/r^), r^ Moliere radius. From two neighbouring points x. , x., 
the lateral age parameter s^^ in given by ^3 


where 


s. . = Ln (F. .X^.Y^:^)/Ln (X. .Y. .) 
ij ij ij ij ij ij 


■'id ■ . Xjj = r./r. , 


Y, . = (x. +1)/ (x . + 1) . 
^ J i J 


If r. r, , s. . s(r) when r = (r.+r.)/ 2 . 

1 J X 3 

Different behaviour of lateral structure function were proposed 
from different analytical treatment of diffusion equations ( 3 ) or 
Monte-Carlo simulation of e.m. cascades (4). A convenient formula 
was advanced consisting to use NKG formula with a Moliere radius 
reduced by a factor m = O. 78 -O. 2 I s^ (for individual e.m. cascade). 
We observed for single e.m. cascade^that s(r) had, versus (3) or 
(4) a behaviour similar to fig. 1 ; incorporating those results in 
EAS- 3 D simulation, we ascertained that this typical behaviour of 
s(r) survived in EAS lateral distribution for all sizes and levels. 
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Recently we developed a very detailed Monte-Carlo 3D-simulat ion , 
including all possible causes of deviation (for instance, the dou- 
ble body decay of each is completely described to produce 

energy and director cosines of each outgoing y, < is correla- 

ted with central rapidity density...) ; the model used for nuclear 
interaction is the multicluster phenomenological one described in 
HE 4.1-9,10. A surprising agreement is obtained with Akeno data (5) 
in favour of the parametrizat ion (fig. 1) 

2 

s(r) = ai + + y , where ^ = Ln (r/r^) 

For a uniform density of detectors, the answer of an experimen- 
tal array will be averaged on all the shower disk seen by the array. 
The radius R of this shower disk depends on the detector area (or 
density threshold) and it comes : 

i = ^ J s(r) dr = (aF^(X) + j3G^(X)) + y with X = R/r^ 

where F^(X) = XLn^ X - 2X(Ln x -1) , G (X) = XLn X - X. 

According to' the shape of fig* 1, any attempt to estimate s in 
[o - 50 m] and [50 - 150 m] will give the same value of the inte- 
gral (6). 


3> Comparison with experiment 

The variation of s versus size has been calculated for MPM 
(nucleon primaries, 6 = 22°5) at Akeno level 2 (fig. 2). Four values 
of R have been taken corresponding to 4 e /m , 2 , 1, 0.5 e~/m . The 
first value corresponding to 0.25 m^ detector is in very good 
agreement with the experimental data. At Tian-Shan level, we have 
plotted the correlation between s^ and s^ , the first parameter 
being estimated by NKG function from numerical values simulated at 
6 and 70 m (fig* 3 ) for different sizes between 10 ^ - 10 ^ particles. 
A general correlation appears : s^ = 1.434 " 0.243. 

The correspondance giving s at Akeno and Moscow altitu^cjes has 
been also obtained as s^ = 1.246 + 0.044 and s^ = 1 . 157 s^^+ 0 . 183 . 

We note that with the present assumption included for e.m. cascade, 
the correlation is not independent on level (fig. 4). 


4. Discussion 

The behaviour versus r of s(r) is not smeared out by the hadro- 
nic cascade and survives in EAS (fig. 1). The agreement obtained 
previously with high multiplicity model is now obtained with the 
multicluster phenomenological model ^ 

(describing p-p data) up to 5*10 GeV (limit of our Monte- 
Carlo simulation) as well at Akeno for s and s(r) as in Tian-S^an 
for s^ _ . The age parameter data doesn’t suggest between 2.10 
5 . 10 ^ GeV any increase of primary mass and supports better a 
nucleon dominance in primary cosmic rays. 
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